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ABSTRACT 

Tidal disruption of main sequence stars by black holes has generally been 
thought to lead to a signal dominated by UV emission. If, however, the black 
hole spins rapidly and the poloidal magnetic field intensity on the black hole 
horizon is comparable to the inner accretion disk pressure, a powerful jet may 
form whose luminosity can easily exceed the thermal UV luminosity. When the 
jet beam points at Earth, its non-thermal luminosity can dominate the emit- 
ted spectrum. The thermal and non-thermal components decay differently with 
time. In particular, the thermal emission should remain roughly constant for 
a significant time after the period of maximum accretion, beginning to dimin- 
ish only after a delay, whereas after the peak accretion rate, the non-thermal 
jet emission decays, but then reaches a plateau. Both transitions are tied to a 
characteristic timescale tEdd at which the accretion rate falls below Eddington. 
Making use of this timescale in a new parameter-inference formalism for tidal 
disruption events with significant emission from a jet, we analyze the recent flare 
source Swift J2058. It is consistent with an event in which a main sequence 
solar-type star is disrupted by a black hole of mass ~ 4 x IO^Mq. The beginning 
of the fiat phase in the non-thermal emission from this source can possibly be 
seen in the late-time lightcurve. Optical photometry over the first ^ 40 d of this 
flare is also consistent with this picture, but is only weakly constraining because 
the bolometric correction is very uncertain. We suggest that future searches for 
main sequence tidal disruptions use methods sensitive to jet radiation as well as 
to thermal UV radiation. 

Subject headings: accretion, black holes 
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Introduction 



Now th at we know that supermassive 



of galaxies ( iKormendy fc Richstone 



1995 



jlack holes frequently can be found in the centers 



Giiltekin et al. 



20091 ) . we can also expect that 



ordinary stars occasionally pass near enough to them to be tidally disrupted. Theoretical 
estimates of the rate a t which this might happen range fror n 3 x 10~^Mp|f^ galax y"^ yr~^ 



( iBrockamp et al. 



20111) to 6.5 X 10-^M-]^f galax y-^ vr"^ flWang &: Merritt 



observational estimates are in the same bal 



Maksym et al 



2010|), although 



Donley et al 



p ark (Gezari et al. 



2008 



2004) 



Esquej et al 



. Most 



2008 



( 120021 ) estimated a rate of X-ray flares as low 



as ~ 10 ^ galaxy ^ yr ^. 



More than twenty years ago, 



ReesI (Il988[ ) outlined the basic dynamical processes 



relevant to these events. The most likely trajectory for a tidal disruption is one in which 
the star approaches on a nearly parabolic orbit. When the star is disrupted, its different 
parts acquire binding energies that span a very large range. The material of the star begins 
returning to the vicinity of the black hole after a delay corresponding to the orbital period 
of the most tightly bound matter. If the disrupted star has a uniform distribution in orbital 
binding energy per unit mass across that large range, matter continues to fall in at a rate 
dM/dt oc ( Phinneyl 119891 ). This matter is captured into an accretion disk. Within this 
disk, gravitational energy is dissipated and the heat is radiated in the usual qua si-therrn al 



fashion. At the peak accretion rate, the luminosity would be super-Eddington (jUlmer 
19991 ): if it is thermally radiated by a disk of area corresponding to a few gravitational radii 
of the black hole, the associated temperature would be in the EUV or perhaps the soft 
X-ray band. Consequently, most searches for such events hitherto have been carried out in 
the EUV region. 



Many o 



' these expectations were u pended by the remarkable object Swift J1644-I-57 



( iLevan et al. 



2011 



Burrows et al 



201 ll ). During its time of peak luminosity (the first ~ 2 d 
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of the outburst), its lightcurve, with its sequence of flares lasting ~ 1000 s alternating 
with quiescent periods ~ 5 x 10*^ s long, bore little resemblance to the prediction of a 
rise-time of ~ 10® s followed by a smooth roll-over to a t~^/^ fall-off. Moreover, whereas all 
expectations were for the emission to peak in the EUV, the greatest part of the luminosity 
was carried in photons with energies ~ 100 keV, indicating electron energies far greater 
than the expected ~ 100 eV, while the shape of that spectrum {vFy gradually rising from 
2 keV to 100 keV) certainly did not suggest any variety of thermal emission. In fact, the 
most successful emission models i nvoke non-thermal s ynchro- Compton radiation from a 



relativistic jet (IBloom et al 



2011 



Burrows et al. 



201lh . 



If most of the radiation actually comes from a relativistic jet, then clearly some of 
the conventional assumptions must be revised. The observed radiation has nothing to do 
with quiescent radiative cooling of an optically thick accretion flow. Indeed, if the jet 
is powered by something resembling the Blandford-Znajek mechanism, its real source of 
energy is rotation of the black hole; the only function of the mass accretion is to support 
the magnetic field threading the black hole's horizon. 



We have recently shown (IKrolik fc Piran 



20 111 ) that a possible solution to all these 



conundrums is that this event was caused by tidal disruption not of a main sequence star, 
but of a white dwarf, and the disruption takes place over several passes, not in a single go. 
Most importantly for the question examined in that paper, a white dwarf disruption results 
in depositing so much mass in such a small region that the accretion flow is immensely 
optically thick. As a result, the matter's cooling time is far longer than its inflow time. 
In such a state of affairs, thermal radiation from the accretion flow is suppressed, and the 
disk maintains a high pressure, one capable of confining a strong magnetic field close to the 
black hole. 



Our goal in this paper is to inquire whether a similar state of affairs could occur when 
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the victim of tidal disruption is a main sequence star. We estimate here both the thermal 
emission of the disk and the jet power. The first is estimated using standard accretion disk 
theory. The latter is estimated in terms of the Blandford-Znajek mechanism, using the 
disk pressure near the ISCO as a measure of the black hole's magnetic field. Although a 
tidal disruption of a main sequence star takes place on a much larger lengthscale than the 
disruption of a white dwarf, for most of the relevant parameter space a similar situation 
holds: the accretion is super-Eddington, thermal radiation is suppressed, and conditions for 
emergence of a strong jet are established. 



Our work is complementary to two related efforts. 



Giannios fc Metzeerl (1201 If ) 



investigated possible radio emission from a jet receiving a fixed (small) fract i on of the 



Lei fc Zhand (120111 ) have 



accretion energy released by accreting tidally disrupted matter, 
suggested a similar picture, but approach it rather differently. In particular, they use thin 
disk approximations for both the sub- and super-Eddington regimes, their scaling with 
black hole mass does not include the relation between disk thickness and accretion rate in 
the radiation-dominated sub-Eddington phase, and they do not discuss the luminosity of 
the thermal disk. We begin in § [2] with a brief discussion of tidal disruption physics. We 
then discuss accretion dynamics in this context and the jet and disk outputs in § [31 In 
§ m we show how this approach can be used to constrain a number of otherwise- unknown 



parameters of tidal disruptions and apply this method to Swift J2058 (jCenko et al 



2Qll|), 



a second example of a jet-dominated tidal disruption. We summarize our results in § [51 



2. Tidal Disruption 



In order for the tidal gravity of a black hole to strongly affect a star, the pericenter of 
the star's orbit cannot be much larger than 



Rt ^ m/ffl^MT-^M-lilR,, 



(1) 
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where k is the apsidal motion constant (determi ned by the sta r's radial density profile) 



and / is its binding energy in units of GMl/R^ (jPhinney 



1989I ). Here A^^, is the mass of 



the star in solar units, Mbh,6 is the black hole's mass in units of IO^Mq, and Rg is the 
black hole's gravitational lengthscale GMbh / ■ We have approximated the main sequence 
mass-radius relati on by -R^ ~ Rp,M^ £ ^ 0.2 for 0.1 < A^* < 1, but increases to ~ 0.4 



for 1 < < 10 (iKippenhahn fc Weigert 



19941 ). The analogous expression for non-main 



sequence stars can be found by substituting the appropriate mass-radius relation; for 
example, for white dwarfs with masses well below the Chandrasekhar mass, ^ = 4/3, and 
the radius of a IMq star is ^ 0.011 -Rf^. The ratio k/ f c:^ 0.02 for radiative stars, but is ^ 0.3 



for convective stars jPhinney 



19891 ). Although the ratio k/f appears in the expression for 
Rt raised to only the 1/6 power, and would therefore appear to be an innocuous correction 
factor of order unity, it can be quantitatively significant because the timescale for the flare 
is oc R^ and the contrast between the values oi k/f for radiative and convective stars is 
^ 15. In numerical estimates, we scale k/f to the value for fully radiative stars because 



this is a reasonable approxirn ation for main sequence stars with O.4M0 < M < 10M( 



( iKippenhahn &: Weigert 



19941 ). 







The work done by the tides stretching and rotating the star is removed from the 
orbit, so on average, the star's material becomes bound to the black hole with a specific 
binding energy of order the star's original specific binding energy ~ GM^/R^. Because the 
semi-major axis a of an orbit with binding energy Eb is GMbh/{'2Eb) when <C Mbh, 
the semi-major axis of material with the mean energy is much larger than Rt, 



a^ean ~ (1 / 2) {k / f)-'/%MBH / M^f^' Rt- 



(2) 



However, when the star is disrupted, its different parts acquire binding energies that span 
a very large range. Suppose that the actual pericenter is Rp = PRt, where /3 ^ 1 is called 
the "penetration factor" . The most bound matter may have a binding energy as great as 
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GMbhR*I Rp, leading to a semi-major axis only 



(3) 



much smaller than amcan- Tidal torques during the initial stage of the en counter cou ld also 



19881); in 



lead to a somewhat larger binding energy for the most bound material (iReesI 
terms of the discussion in this paper, the effects are indistinguishable from those due to a 
decrease in {3. 

On the assumption that the disrupted star has a uniform distribution in orbital binding 
energy per unit mas s, matter retur ns to the region ~ Rp from the black hole at a rate 



dM/dt oc {t/to)-^/^ flPhinnev 



19891 ). The characteristic timescale to for initiation of this 



power-law accretion rate is the orbital period for the most bound matter 



-Porb ( '3'min 



1/2 



(4) 



the period fo r an o rbit with semi-major axis amm- More detailed calculations (e.g.. 



Lodato et al. 



(120091 )) indicate that, almost independent of the star's internal structure, the 
mass return rate does follow this behavior at late times, but at early times the rate rises 
rapidly after a delay ~ -Porb(amin) and then gradually rolls over to the t~^^^ proportionality, 
transitioning more slowly when the star's pre-disruption structure is more centrally 
concentrated. 

In this picture, once the matter returns to the vicinity of Rp it is captured into an 
accretion disk whose inflow time tin ^ -Porb(a.mm)- As the matter moves inward through 
this disk, there is local dissipation of the conventional accretion disk variety, and the heat 
is radiated in the u sual quasi-thermal fashion. At the peak accretion rate estimated by 



Lodato et al. 



(120091 ) the luminosity would be 



-1/2 



(3 ^ erg/s 



(5) 
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for radiative efficiency 77. In Eddington units, this luminosity is ^ 800/3 ^ ( r; /0 . 1 ITW 1^"*"^^^ ^Mp|/^ 



for radiative stars and ~ 4 times smaller for convective stars (see also 



Strubbe fc Quataert 



Ulmeil (|l999|); 



(120091 )): as we will discuss in § [HI super-Eddington accretion may 



result in only a fraction of this inflow rate reaching the black hole. If such a luminosity were 



thermally radiated by a disk of area ~ O{10)nR'^, the associated temperature would be 



char 



6 X 10^M^}+^^^/^M, 



-5/8 



k/f 



13-'/^ K, 



implying that most of the radiation would emerge in the soft X-ray band. 



(6) 



3. Accretion dynamics and the jet and disk powers 



We begin by estimating conditions at the time of peak accretion rate. iLodato et al. 



( I2OO9I ) estimate that ~ 1/3 of the star's mass should arrive at Rp after a delay ~ -Porblamm)- 
Because M falls oc t~^^^ thereafter, dM/dlnt is greatest during this period, so this is both 
the period of greatest luminosity and the period in which most of the energy of the entire 
event is radiated. As the matter streams in, its orbital velocity is roughly the free-fall 
velocity at Rp. Collision with any mass in a more circular orbit involves velocity differences 
of order the orbital speed, so the post-shock temperature is comparable to the virial 
temperature. Following conventional accretion disk theory in the supposition that the inflow 
rate is controlled by angular momentum transport and writing the vertically-integrated 
stress as some number a times the similarly vertically-integrated pressure, we estimate the 



inflow time at Rp as 



ti„ ~ a-\Rp/H)\RT/RGf'\R,/c) {^^"' P'" W 
~ 2 X l{}\a/Q.lY\k/fY/\Rp/HfMl-'^^'^~^^0^'^ s. (7b) 

The quantity H is the disk thickness, which cooling may reduce to be less than i?p, while 
a is the usual ratio of stress to pressure. Note that tin depends only on M*, and not at 
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all on Mbh, because it is fundamentally a dynamical time, and at the tidal radius that is 
determined by the dynamical time of the star ~ {Rl/GM^y^"^. Because tin ^ -Porblc^mm, 
the mass held at Rp is determined by a balance between the slowly-changing rate at which 
returning tidal streams deliver mass and the more rapid inflow rate. The characteristic 
Thomson optical depth from the middle of the material to the outside at Rp is then 

-1/12 



~ 3 X 10\a/0.1)-\Rp/HYM'^^^''^^/'M-Jil ( ^ 



0.02 



/3- 



-7/2 



and the cooling time is 



''cool 



TrH/c ~ 4 X 10^(«/0.1)-i(i?,/i7)A^^6«)/^M-^/^ ( 



0.02 



(9) 



The effectiveness of cooling can be measured by comparing tcooi to the inflow time tj 
at that radius: 

-1/6 



tcooi/tin ^ 500{H/Rp)M;'/'^^'/'^^M-'Jl 



k/l 
0.02 



/3" 



(10) 



If the initial temperature of the matter arriving in the disk is close to the virial temperature, 
it cools efficiently only if the black hole mass is ~ 10^ Mq or the star's mass is ~ O.IMq. 
Although it is possible for stars on specific trajectories to be tidally disrupted by spinning 
black holes with masses up to ~ IO^Mq, the probability of tidal disru ption (as oppo sed to 



direct capture) declines rapidly when Mbh,6 ~ lOOA^i '■^^^^^ or more (iKesdenI 
this reason, unless the matter arrives at Rp substantially cooler, we expect that at the time 
of peak accretion rate, the inflow will never be able to cool efficiently. 

Another way of viewing this result is to n ote that photon trapping is expecte d 



201 ih . For 



whenever the accretion rate is super-Eddington (IBegelman 



1979 



Abramowicz et al 



19881 ). 



As we have already seen, the characteristic peak luminosity in these events exceeds the 



Eddington luminosity by a factor ~ 800(77/0.1)A^' 



(1+30/2^-3/2 1/2 



BH,6 I 0.02 



/3- 



Thus, to 



order-of-magnitude accuracy, the dividing line between sub- and super-Eddington accretion 



-lo- 



in main sequence tidal disruption events falls at roughly the same place as the dividing 
line between the radiatively efficient and photon-trapping regimes, and both are near the 



maximum b 



( lUlmer 



1999 



ack hole mass at which the re is a significant probability of tidal disruption 



Strubbe fc Quataert 



20091). 



3.1. Jet power 



Rotating bl ack holes whose horizons a r e threaded by large-scale po l oidal field can driy e 



relativistic jets ( iBlandford fc Znajek 



1977 



McKinney fc Gammie 



2004 



Hawley &: Krolik 



20061 ). At the order of magnitude level, the Poynting luminosity of such a jet can be 



estimated through a simple dimensional argument: it must be ~ cB^Rg, where is 
the poloidal field intensity on the horizon. The actual magnitude of the luminosity can 
then be written as Lj^t = f{ci'/M)c{B'^/87r)R'^, where f(a/M ) is a dimension l ess fu nction 



Beckwith et al. 



fl2009[ ). who 



that should increase with \a/M\. To estimate B^, we follow 
demonstrated that the magnetic pressure near the horizon is generally bounded above by 
the midplane total pressure near the ISCO and bounded below by the magnetic pressure at 
that location. Our first task is therefore to estimate the disk pressure near the ISCO. 

When the black hole mass is very large (i.e., roughly the maximum black hole mass 
capable of disrupting the star before it enters the black hole), the infiow can be both 
sub-Eddington and cool efficiently even at the time of peak accretion rate. If tcooi/^m < 1 at 
Rp (eqn. [TOj) . the matter can cool to the temperature at Rp associated with steady accretion 
at the prevailing rate. At smaller radii, the temperature profile of the accretion fiow, as 
well as its o t her pr operties, can be described by standard stationary disk solutions (e.g.. 



as m 



KrolikI (119991 )). Although the accretion rate is sub-Eddington, it is nonetheless high 



enough that the fiow will certainly be radiation-dominated. The pressure in the midplane 
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of the inner disk is then 



2^ 



Pmid 



anR, 



iR/Rg)-^/^X{R/R 



■3/' 



(11) 



where k and kt are the actual and the Thomson, opacity, respectively, and the function 
X{R/Rg) is a place-holder to account for how these scalings are altered as the ISCO is 
approached. Note that pmid in the sub-Eddington radiation-dominated regime is inversely 
proportional to Mbh- Although the outgoing thermal radiation flux is proportional to the 
accretion rate, the optical depth in this regime is inversely proportio nal: as a result, the 



midp lane radiation pressure pmid is independent of the accretion rate (IModerski fc Sikora 

mm- 



Conversely, in the low mass limit (which applies to the majority of relevant black holes), 
the inflow is both super-Eddington and radiatively inefficient because photon-trapping 
prevents effective radiative cooling. In this regime, the effects of strong radiation forces and 
possible radiation- driven outflows make the t heory less clear ( consider, e.g., the range of 



views p resented in 



fl2009h 



Abramowicz et al 



Dotan fc Shavivl (120 llh 



19881): 



King fc Poundsl (120031 ): 



Strubbe fc Quataert 



Begelmanl (1201 ll )). For our purposes, it is sufficient to 



estimate the pressure using dimensional analysis and simple scaling arguments. Precisely 
because the gas cannot cool itself, the ratio of total (i.e., radiation plus gas) pressure to 
density remains at a level close to the virial temperature as the flow moves inward, so that, 
unlike the sub-Eddington radiation pressure-dominated regime, the pressure when the flow 
is super-Eddington is proportional to the accretion rate and the flow conflguration is nearly 
round. The midplane pressure (again, dominated by radiation) is then 



Pmid 



(12) 



where Cg is the effective sound speed (including radiation pressure) and Worb is the speed 
of a circular orbit at that radius. Here we have substituted a for the ratio of inflow speed 
v^. to orbital speed because fr/^Vb ~ <^ when the disk is geometrically thick. We deflne fn 
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to be the accretion rate delivered to the outer edge of the disk normahzed to Eddington 



makes the radiative efficiency difficult to estimate. When the distinction is important, we 
will use mo for the value of this normalized accretion rate at its peak value. The parameter 
q gives the fraction of m arriving at the black hole; when m > 1, it is possible for g < 1, but 
when m < 1, in general we expect = 1. Modulo corrections of order unity, when electron 
scattering dominates the opacity and the relevant radii are ~ -Rg, the only difference 
between this expression and the previous one is that this one is also oc qm. 

Because our estimates for the total pressure in both the sub- and super-Eddington 
regimes are so similar, it is convenient to write an estimate for the jet luminosity in a single 
form: 



where a should be interpreted as the usual stress parameter in the sub-Eddington regime, 
but only as tv/^orb in the super-Eddington case. When MBH,Qf{.o./M)/{al3h) = 1 and (in 
the super-Eddington regime qm = 1), the magnitude of the jet power is 1 x 10^'^ erg s~^. 
The quantity Ph is the ratio of the midplane total pressure near the ISCO to the magnetic 
pressure in the black hole's stretched horizon. In very rough terms, we might expect 
PhC( ~ 0.1-1. Thus, for fixed black hole spin, the power in the jet is a constant fraction of 
the Eddington luminosity, independent of accretion rate, in the sub-Eddington regime, but 
rises in proportion to the accretion rate in the super-Eddington regime. The reason is that 
at lower accretion rates, the radiative losses curb any increase in pressure with increasing 
accretion rate, whereas once the accretion rate becomes super-Eddington, the photons are 
trapped in the accretion flow. 

However, there are several subtleties hidden in this estimate. The ffist is that f{a/M) 



units with unit efficiency, i.e., m = Mc^/Le because the uncertain effect of photon-trapping 




(13) 
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is not well-determined. iBlandford fc Znajekl (119771 ) worked out an expa n sion i n small a/M 



for time-steady split-monopolar field distributions. 



Tchekhovskoy et al 



teoioh developed a 



similar expansion in terms of the black hole rotation rate. Both required t he fie. 



Tchekhovskoy et al. 



d rotation 



mm found 



rate to be exactly half the black hole rotation rate, 
that their expression compared well to axisymmetric force-free MHD simulations with an 
imposed magnetic flux function. Unfortunately, neither axisymmetry nor an imposed flux 
function may be a good approximation to the behavior of time-dependent 3-d relativistic 
jets whose field configuration is created self -consistently by a.ccreti on dynamics, which, in 



turn, may be influenced by black hole spin (IHawley fc Krolik 



2006f ). It is therefore difficult 



to be confident about much more than that jet power should increase with faster black hole 
rotation. 

The second subtlety is that the relevant field intensity on the horizon is the time-average 
of the poloidal component. When the sign of the vertical field through the equator ial plane 



20081). Tidal 



changes frequently, this intensity can be strongly suppressed (iBeckwith et al.l 
disruptions, however, may be a particularly favorable case for jet-launching because the 
extreme elongation of tidal streams might create very large-scale structures in the field with 
exactly the sort of coherent direction required to support powerful and long-lived jets. 

Third, up to this point we have been careful to discuss only the two extreme limits: 
high-mass, sub-Eddington, radiatively efficient; and low-mass, super-Eddington, and 
radiatively inefficient. That is because the transition between them is not well understood. 
For estimating both the jet luminosity and (later) the emergent thermal luminosity, the 
key parameter is the ratio tcooi/^in- Several problems must be solved before its transition 
between the two extreme limits can be defined. There may be parameter regimes in 
which this ratio is > 1 at Rp if the gas is deposited there with T ~ Tvir even though the 
temperature in a steady-state flow at the specified accretion rate is sufficiently cooler that 
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the gas could cool in less than an inflow time. On the other hand, there are accretion rates 
for which the flow in steady-state is radiatively efficient at Rp, but not at smaller radii. 
The quantitative character of the transition is also sensitive to exactly how the dissipation 
profile varies with radius near the ISCO. In view of these uncertainties, we will adopt here 
the simplest solution, understanding that it is only provisional: we will define the transition 
point by the intersection between the two expressions for Ljot; we will adopt a similar 
prescription for the thermal luminosity. 

Fourth, there is the question of what fraction of the returning mass actually reaches 
the vicinity of the black hole, i.e., the correct value of q. This is important in the 
super-Eddington regime where Ljet oc qm. If a significant fraction of rri is ejected as a wind 
(as could well happen during the super-Eddington period), q could be rather small, and the 
jet power would be reduced accordingly. 

Lastly, jets in Galactic b 



( IRemillard fc McChntock 



ack hole binaries turn off when m rises toward ~ 1 
20061 ). On the other hand, the existence of radio-loud AGN with 
strong optical/UV continua and emission lines suggests that larger black holes can somehow 
support jets even when the accretion rate is near Eddington. Luminous quasars are 
thought to accrete a t rates such that 0.01 <^ L/ Le ^ 1 (the lower end of this range favored 



by, e.g. 



Kelly et al 



(120101 ) ■ who estimate black hole masses frq r n the 



Liu et al 



quasar luminosity 



(120091), whose black hole 



and broad-line width; the higher end favored by, e.g. 
masses come from the bulge dispersion in obscured quasars). Although these statistics 
are dominated by the radio-quiet variety, they likely apply to radio-loud quasars as well 



because there are on 



radio-quiet quasars (IRichards et al. 



y slight differences b etween the optical/UV continua of radio-loud and 



201ll ). a fact suggesting that the inner disks of these 



two categories of quasar are very similar. A more direct indication comes from Narrow 
Line Seyfert 1 galaxies (NLSls). These objects are thought to accrete at near- or possibly 
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super-Eddington rates (IBoUer et al.lll996l : iBorosonI 12002 



it has been found that a fraction of them are radio-loud fKomossa et al 



Wang fc Netzei 120031) . Recently, 



20061 ). and several 



20091 ): such objects must 



of these radio-loud cases are strong 7-ray sources (lAbdo et al. 
certainly have strong jets. Because the black holes in galactic nuclei responsible for tidally 
disrupting stars have masses closer to those in AGN than those in Galactic black hole 
binaries, we will assume that the AGN example is the guiding precedent here. 

In spite of all those caveats, we believe that the fundamental characteristics of this 
system make production of a significant relativistic jet very likely when the black hole 
spins rapidly. Moreover, if there is any jet associated directly with the black hole, on 
dimensional grounds its luminosity must scale with cpmid^g, although the dimensionless 
factor multiplying this quantity could be far from unity. 

The previous estimate for jet power (eqn. [T3|) posed the issue in terms of the accretion 
rate relative to Eddington. However, this accretion rate is predicted (approximately) by the 
dynamics of tidal disruption: 



Ljet~lxlO^ — 3— erg/s ^ ^ ^ 



BH,6 



where 



M 



BH let 



4 X lQ\m/mof"'M^}+^^'^ 



k/l 
0.02 



-1/3 



Mbh Mbh jet, 
(14) 

(15) 



In other words, the jet luminosity has a minimum as a function of Mbh- When the black 

— 1/2 

hole mass is smaller, Ljet oc M^^ because that is the photon-trapping regime, making 
p oc gm, and tidal disruption mechanics make rriQ oc M^'^'^; when the black hole mass is 
larger, the accretion flow is in the conventional sub-Eddington regime, and Ljet Mbh 
because p is independent of M but oc M^]^, while the black hole area oc M^jj. The 
transition mass between the two regimes, MBHjct, decreases with decreasing accretion 
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rate. Thus, a system that is initially in the photon-trapping regime will eventually become 
sub-Eddington later. On the other hand, a large black hole system that is sub-Eddington 
initially will remain so throughout the event. Even if only a small fraction of the fallback 
accretion rate reaches the black hole (i.e., g ^ 1), the jet luminosity would still increase 
toward smaller black hole masses provided q scales less rapidly with black hole mass than 

It is important to recall that at the time of peak accretion, the black hole mass dividing 
the low- and high-mass regimes is almost as large as the yery la rgest black hole mass 



permitting tidal disruption of main sequence stars (iKesdenl 1201 ll ). This conclusion would 
be weakened somewhat for stars with large internal convection zones (stars at either the 
low or high end of the mass distribution), but MsHjct decreases by only a factor ^ 2.5. 
At the peak of the flare, therefore, we expect essentially all events to be effectively in the 
low-mass, super-Eddington regime. Only at later times, as the accretion rate falls, does the 
high-mass regime become relevant. 

Note also that the jet luminosity is not directly limited by the Eddington luminosity 
because its energy is not directly related to the accretion flow. The pressure in the accretion 
flow confines the magnetic field, but the magnetic field taps the rotational kinetic energy 
(i.e., the difference between the total mass and the irreducible mass) of the black hole. For 
the same reason, the jet efficiency measured in rest-mass units, which is f{a/M)/{47rPhCt), 
is likely to be ~ O(O.l), but can in principle be greater than unity. 

It is also important to emphasize that the jet power does not translate directly into the 
photon luminosity we see. On the one hand, the photon luminosity of the jet is generally 



a small fraction ?7jct of the initial P oynting luminosity. 



jets, for example, ?7jct ~ 10 ^-10 ^ (jCelotti fc Ghisellini 



n the conditions relevant to blazar 



20081 ). On the other hand, whatever 



photon luminosity emerges will in general be strongly beamed due to jet coUimation and 
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relativistic kinematics: B ~ min(47r/Af2, 2r^), where AQ is the sohd angle occupied by the 
jet and F is its Lorentz factor. The beaming factor B can easily be ~ 10^. 



3.2. Thermal luminosity 

The thermal photon luminosity behaves differently as a function of black hole mass. It 
can be written as //therm = VI'^^^e, where we expect that t] depends on spin for m < 1, but 
declines once m is large enough that the portion of the accretion flow where the majority of 
the luminosity would otherwise emerge enters the photon-trapping regime. In this limit, the 
photon luminosity is capped at ^ Le- As we have already discussed in the context of the 
jet power, the critical m above which this occurs is not well-determined. Thus, we estimate 



/^therm =i 1-5x10 erg/s 



Mbh,6 Mbh ^ MBHthcrm 

800(r//0.1)(m/mo)A^l'+'^^/'M-i/^ (^^j"'^'/3-3 Mbh A/Bmhcrm- 



(16) 

In the sub-Eddington regime (the case of larger black hole masses), the peak emerging 
luminosity is exactly the characteristic luminosity we estimated previously (eqn. |5]). The 
black hole mass dividing the super- and sub-Eddington regimes is 

MBm^err. = 9.0 X 10^(V0.1)^/3(m/mo)^/^Ml^+^^)/^ '^^-^M^. (17) 

In other words, Ltherm has a maximum at a black hole mass ^ MBHthcrm- The thermal 
luminosity is oc Mbh at smaller masses because photon trapping limits the emergent 

— 1/2 

luminosity to L^;; it is oc M^^j at higher masses as it reflects M in the radiatively efficient 
regime. Much as in the case of the jet luminosity, the low-mass case will almost always 
be the regime relevant to the peak of the flare, but the high-mass case may be seen as 
rh decreases. Note that our critical mass for phot on-trapping, MBHtherm, is several times 



greater than the mass estimated by 
factor in the expression for Rj-. 



Ulmerl (119991 ) because we have allowed for the k/f 
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These trends for the peak values of Ljct and Ltherm as functions of Mbh are illustrated 
in Figure [U The Ljet in this figure is to some degree a maximal estimate, as it supposes 
both q = 1 (i.e., little loss of accretion rate due to outflows) and f{a/M) = 1. However, 
as the figure shows, if the black hole rotates rapidly enough to make f{a/M) not too 
small, for Mbh ^ lO^M©, a very large proportion of the returning mass would have to be 
ejected in order for Lj^t to fall to a level at which it is merely comparable with Ltherm- For 
rapidly-spinning black holes, Lthcrm ~ Ljct only when the mass is the largest permitting any 
tidal disruptions at all. 

Thus, the peak output is almost always dominated by the jet unless the black hole 
rotates rather slowly and the black hole mass is ~ lO^M©, a mass so large that the tidal 
disruption probability is considerably less than unity. However, whether non-thermal 
radiation from the jet or thermal radiation from the disk dominates the observed spectrum 
depends on the interplay of the jet radiative efficiency and beaming with the intrinsic 
jet/disk ratio determined by a/M and Mbh- The jet is favored by viewing angles along its 
axis, low black hole mass, high spin, and large ?7jct; the disk is favored by viewing angles 
outside the favored cone, high black hole mass, slow spin, and small ?7jet. 



3.3. Time-dependence 



Now consider what happens at later times, as matter of progressively smaller binding 
energy returns to the disk. If the mass per unit energy of tidally disrupted matter 



is independent of binding energy, t he accretion rate into the disk at radii y Rp 



oc (t/to) where to ~ -Porb(amin) (iRees 



1988 



Phinney 



19891) 



Lodato et al. 



alls 



(120091) found 



that, depending on the degree of density concentration in the star, the mass accretion rate 
might fall more shallowly than this when t ^ to, but gradually tends toward t~^^^ when 
t ^ tf). If nearly all the matter returning to the disk makes it all the way to the black hole. 
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Fig. 1. — Peak jet power (blue) and thermal (red) luminosity as a function of black hole 
mass, following the prescriptions of equations [H] and [lEl All scaling factors except Mbh,6 
in those expressions are set to unity. 
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the time-dependence of the accretion relevant to the disk and jet powers should match the 
time-dependence of the accretion rate reaching the outer disk; if, however, super-Eddington 
conditions lead to a significant fraction being expelled, the accretion rate reaching the event 
horizon falls more slowly because the expelled fraction can also be expected to diminish. 

When m first begins to decrease, the thermal radiation hardly changes because 
photon-trapping continues to limit Lthcrm to Le- Consequently, in the period near and 
shortly after the peak luminosity, the thermal light curve should be almost flat and only 
slowly roll over toward {t/to)"^/^. 

On the other hand, Ljet in those sources is oc gm, so the jet power falls, provided the 
accretion rate reaching the black hole is truly super-Eddington. The rate of decrease may 
be slower than if q increases as m falls. With declining m, the black hole mass at 
which the flow switches from non-radiative to radiative also decreases. Consequently, even 
though virtually every source begins in a non-radiative state, those in which the black hole 
mass is relatively large may ultimately become efficient radiators at later times. When that 
changeover occurs, the jet power remains constant, while the thermal luminosity begins 
to fall. This timescale, when the accretion rate and therefore the thermal luminosity fall 
below Eddington, we call tEdd- It should be emphasized that all of these remarks pertain 
to the expectation value for the jet luminosity; relativistic jets are generically unsteady, so 
fluctuations at the order-unity level should be expected around all these trends. 

Figure [2] gives a schematic view, beginning at the time of peak accretion rate, of what 
might be expected in terms of the light curves for the jet power (before allowance for 
beaming and radiative efficiency) and the thermal luminosity. For the parameter values 
chosen {Mbh^ = 10, all other scaling parameters unity), Ljet falls to the level of Ltherm at 
almost the same time, t ^ 7to, as Ltherm enters the sub-Eddington regime and also begins 
to decline. From that time to t ~ 30to, both fall together, maintaining similar power levels. 
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Finally, after t ~ 30to (i-e., a time larger by r] 3/5^4 than the time at which the thermal 
luminosity begins to decline), the jet luminosity stabilizes, while Lthcrm continues to fall. 

The ratio Ljct/ Lthcvm as a function of time for fixed black hole mass can be described 
much more simply. In rough terms, it is just 

-^jet ^ max(l, gm)/(a/M) 
Ltherm min[l , ?7(a/M)m] 

where ri{a/M) is the intrinsic (i.e., without photon trapping) radiative efficiency for a given 
spin parameter. At high accretion rate, the ratio is oc qrh; in the transition range where 
photon trapping is marginal, the ratio may change relatively slowly; at low accretion rate, 
the ratio is oc m~^. In any particular event, m rises rapidly (on a timescale ~ to) at first, 
but then declines slowly. The ratio of observed ]et luminosity to thermal luminosity follows 
the same trend, modulated by any evolution of qBrj-^ct- 

Eventually, the accretion rate falls so low that the disk is no longer radiation-dominated, 
and our estimate for the pressure in the inner disk is no longer valid. Assuming that this 
occurs when the gas pressure becomes comparable to the radiation pressure at r ~ lOr^, the 
critical accretion rate is rh ~ 0.05(r//0.1)-i(a/0.1)"^/SM^]/g. If the accretion rate declines 
as (t/to)~^^^, it happens at a time 

tg ^ 1300to/3-'/'(r//0.1)3/^(a/0.1)3/^°M-|/'°A^f +3«)/io. (^g) 

After this time, although the thermal luminosity continues to decline in proportion to the 
accretion rate, the midplane pressure (and therefore jet luminosity) declines oc m^/^. This 
break time is quite late in the development of the fiare for low-mass events, but when the 
black hole mass is relatively large, it could take place as early as ^ 30to- 



-22- 




Fig. 2. — Jet power (blue) and thermal (red) luminosity as functions of time for Mbh = 
1 X lO^M©. The Eddington timescale tEdd for these parameters is 7to. As in Figured! all 
scaling factors except Mbh,6 and m/rho are set to unity. 
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4. Constraining tidal disruption parameters: the case of J2058. 4+0516 

Although the theory of jet and disk emission from tidal disruptions as we have 
presented it contains a large number of free parameters, in any given event there are also 
potentially a sizable number of observable properties that may be used to constrain these 
parameters. In nearly every event, it is possible to measure the characteristic time to- 
Optical and/or ultraviolet observations often give an indicator of the peak thermal disk 
luminosity, L^o'-, we describe this as only an "indicator" because, as we discuss below, there 
is likely to be a sizable — and uncertain — bolometric correction. When there is hard X-ray 
emission, the peak jet luminosity Ljo can also be obtained; if an event can be followed long 
enough, it may also be possible to measure two times related to the transition from super- 
to sub-Eddington accretion: the time tjet at which the jet luminosity flattens, and the time 
^disk at which the disk luminosity begins to diminish. In this section, we will lay out a 
general formalism for using these observables as parameter constraints and then apply that 
method to a speciflc example. Swift J2058. 4+0516. 



As shown by 



Lodato et al. 



(120091 ). the accretion rate peaks at a time -Porb(amin) past 
pericenter passage (see eqn. [3]), and diminishes thereafter as a power-law in time. If 
we identify the time of peak flare luminosity to with that orbital period, we obtain the 
constraint 

to,.-5.8^^^^^)/x/?,e(^)''%^ (20) 

where to is measured in days. 

Because the peak accretion rate is likely to be super-Eddington for almost the entire 
range of possible black hole masses, we expect the peak disk luminosity to be close to the 
Eddington luminosity. Unfortunately, however, the characteristic temperature of thermal 
disk radiation at the Eddington luminosity is ~ 1 x lO^M^^^g K, indicating that the 
bulk of the light may emerge in the EUV, where direct measurements are very diflicult. 
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Consequently, a bolometric correction that could well be ~ 0(10) or greater must be 
applied to any observable measure of the disk luminosity. Unfortunately, given our current 
limited understanding of disk spectra even when the accretion rate is sub-Eddington, not 
to mention potential dust e xtinction in the host galaxy or possible reprocessing in a wind 



( IStrubbe fc Quataert 



20111), a sizable uncertainty must be attached to any such correction. 



For this reason, we give higher priority to the use of other observables. 

In our description of the time-dependence of the jet and disk luminosities, the transition 
from super- to sub-Eddington behavior occurs earlier for the disk than for the jet. The 
ratio of the accretion rates at these two timescales is ~ t], the radiative efficiency of 
the disk in the trans-Eddington regime. Consequently, the ratio of these two timescales 
primarily constrains the black hole spin, but not any of the other parameters. For the other 
parameters, no additional information is gained by using both timescales; either one will do. 

Suppose, then, that we choose to use the timescale at which the jet luminosity flattens 
as a function of time, tjet- This timescale is typically a few times tEdd- For generality, let the 
accretion rate scale with time as (t/to)~", where n ~ 5/3 is expected. Then we flnd that 

1 /2 

tjet,d ^ 5.8(8 X 103)Vn_;^[l+l/n+35(l/n-l)]/2^(l-3/n)/2 /jSd-l/n)^ ^21) 

where tjet is likewise scaled in days. Combining this timescale constraint with the one 
based on the characteristic flare timescale, we may solve for the black hole mass and the 
penetration factor: 

Mbh,, ^ 4.6 X 10'MJI,%1^ (22) 

and 

1 /6 

/3 ^ 0.093A^p-V3 ~ (23) 

Lastly, using both of these estimates, the (numerous) parameters governing the jet 



luminosity may be constrained: 



B7]i,tqf{a/M)/{Pha) ~ 0.022LjetA6M;H, 



where Ljef,46 is the peak luminosity of the jet in units of 10^^ erg s ^. 
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10"' r 




10 r 
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Fig. 3. — Long-term Swift XRT light curve in the 0.5-10 keV band for J2058 as of 3 January 
2012 (last data point 7 December 2011), combining WT and PC data (data drawn from 



htt p : / / www. swift . psu . edu/ monitoring | . is Swift counts per second. Triangles show la 
upper bounds. The solid line is the result of a x^-squared fit to the lightcurve for t > 7 d. 



20 111 ) may be an 



The recently-discovered fiare source J2058. 4-1-0516 fiCenko et al. 
example of exactly the sort of tidal disruption event to which this formalism is applicable. 
Its lightcurve is shown in Figure El For ^ 10 d, its fiux stayed nearly constant; for the next 
three months, the fiux declined roughly as a power-law in time. Around t ~ 100 d, the 
decline became shallower; unfortunately, the fiux at that point was already near Swift's 
detection limit, so the lightcurve beyond ~ 120 d is a mix of detections and upper bounds 
whose la limits are not very different from the level of some of the last detections. To 
describe this lightcurve parametrically, we fit the data for t > 7 d to a model of the form 
F{t) = At~^ for t < tf) and F{t) = B for t > t^. The minimum in ^-squared was found for 
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n = 1.9 and tb — 120 d; this best-fit model is shown as the sohd hne in Figure |3l That 
is reduced by a value of tb within the span of the data supports our prediction that the jet 
luminosity should become nearly constant after an initial period of decline. 

To use the procedure just outlined, we take to = 10 d, as that appears to be the point 
at which power-law decline begins. Following the results of our fit to the lightcurve, we set 
tjet = 120 d, uncertain as that identification may be, and n = 1.9 rather than 5/3. The 
difference in the parameter inferences due to the latter choice are probably smaller than the 
intrinsic error in the method. Using these numbers, we find 

MsHfi ^ 40A^, (25) 

and 

/3 ^ 0.6>lp-V3 . (26) 

That is, we estimate that the black hole is in the upper range of masses at which tidal 
disruptions of main sequence stars can occur, and the penetration factor is quite modest. 
The expected mass might be a bit smaller if the stellar mass is less than solar; in that case, 
^ ^ 0.2, and /3 would increase. 

The luminosity in this flare appears to be dominated by hard X-rays (photon energies 
^ 10 keV) and peaked (in nominal isotropic terms) at ~ 3 x 10'^'' erg s~^ (or perhaps a 
factor of a few more when contributions from still harder X-rays are considered). The 
constraint on the jet parameters (eqn. [2^ then becomes 

gi3r/jet/(a/M)/ (/3,a) ~ 6M:\ (27) 

For self-consistency, the jet must be reasonably beamed and efficient, and the black hole 
must be spinning rapidly enough to make f{a/M) not too small. If A^* < 1, the beaming, 
jet efficiency, spin, etc. must be somewhat greater. 
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Although the available optical/UV observations are only very rough guides for further 
inference, they are consistent with these values. Optical photometry sugges ts a steeply 



rising spectrum in the observed frame between the g'-hand and the u-band (jCenko et al. 



201ll ). The M-band flux on its own translates to a luminosity uLiy ^ 1.1 x 10^^ erg s~^ given 
the object's redshift [z = 1.18) and assuming a flat cosmology with Hq = 70 km s~^ Mpc~^. 
If Ltherm IS Eddingtou-limited and isotropic, this luminosity alone requires a black hole mass 
7 X IO^Mq. However, the mass estimated in this fashion is proportional to the bolometric 
correction, and, as we have already noted, this may well be at least an order of magnitude. 
If so, the mass required would be in the range just estimated on the basis of to and tjet- 

Lastly, we c onsider what i nform ation may be gleaned from the time-dependence of the 



optical/UV flux. 



Cenko et al. 



( I2OIII ) report g'-hand optical photometry taken within 10 d 
of the flare's initiation and ^ 10 d, ~ 30 d, and ~ 40 d later. They also report w-band 
magnitudes at ~ 10 d and ~ 40 d. In g', there was no detectable dimming until ^ 30 d after 
the flare, but at ~ 40 d, the latest time reported, the flux had diminished by ~ 0.5 mag, 
~ 40%. The M-band flux drops about 0.5 mag from ^ 10 d to ~ 40 d. Thus, in rough terms, 
the optical/UV luminosity appears to have varied very little over the first month or so, even 
though the X-ray luminosity fell by a factor ~ 8; such behavior is in good agreement with 
our prediction that the disk output should remain steady while the jet power falls during 
the super-Eddington phase. If the drop between 30 d and 40 d marks the disk transition to 
sub-Eddington behavior, the factor ~ 3-4 between that timescale and the possible timescale 
of the jet transition is also consistent with our suggestion that the accretion rate at the 
time of disk transition is a factor ~ 77"^ greater than at the time of jet transition. 
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5. Conclusions 

We have demonstrated that tidal disruption of main sequence stars by rotating black 
holes can lead to the launching of a powerful relativistic jet. Because the peak accretion 
rate in most of these events is likely to be super-Eddington, a regime in which photon 
trapping is significant, the thermal luminosity may be held to no more than roughly the 
Eddington luminosity; that is, the thermal luminosity is suppressed well below what the 
accretion rate would predict on the basis of the usual relativistic radiation efficiency. As a 
result, the thermal lightcurve (in bolometric terms) should have a peak that is rather flatter 
than the curve describing the accretion rate. On the other hand, the jet power suffers no 
such suppression, so it can exceed the thermal luminosity of the accretion flow. 

As the accretion rate declines, the jet power diminishes until the sub-Eddington regime 
is reached; after that point, we expect that its power changes little until the accretion rate is 
so low that the disk is no longer radiation-dominated. By contrast, the thermal luminosity 
can be expected to remain roughly constant at roughly the Eddington luminosity as the 
accretion rate falls, diminishing with the classic t~^^^ scaling only after the flow becomes 
sub-Eddington and radiatively efficient. The time at which the thermal luminosity begins 
to fall we have named tEdd! the time at which the jet lightcurve levels out is a few times 

^Edd- 

It follows that in many main sequence tidal disruptions, the jet luminosity may be 
comparable to or greater than the thermal luminosity. Particularly when the black hole 
rotates rapidly, the non-thermal jet photon luminosity may substantially exceed the thermal 
disk luminosity even after allowing for comparatively inefficient radiation of the jet's kinetic 
power. The degree to which we see this non-thermal luminosity depends, of course, on our 
viewing angle. When it is favorable, the non-thermal emission might be significantly larger 
than the thermal. 
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Particularly when the reach of the survey extends only to low redshifts, beaming also 
enhances the accessible population. Although the fraction of sources we see is reduced by B, 
the luminosity distance to which we can see them (for negligible k-correction) is increased 
by jB^/^, implying an increase in accessible comoving volume of B^^'^/{1 + -^max)^, for Zmax 
the greatest redshift at which the luminosity of interest can be detected (provided z^ax ^ 1, 
so that cosmological corrections are small). Thus, if ,B ~ 100 (as from relativistic beaming 
with a Lorentz factor F ~ 10), beaming results in an increase in the number of detectable 
sources by a factor ~ 0(10) provided the population has redshifts smaller than ~ 1. 

For these reasons, it may therefore be desirable to augment tradition al UV-based 



Gezari et al. 



(120081 )) with other 



methods of searching for stellar tidal disruptions (e.g., 
methods more sensitive to jet radiation, perhaps based on the Swift system of all-sky 
monitoring in hard X-rays coupled to 2-10 keV X-ray follow-up (e.g., using possible 
future m issions such as th e Space-based multi-band astronomical Variable O biects Mon itor 



[SVOM 



Paul et al 



Ray et al. 



(I2OIII )] or the Advanced X-ray Timing Array [AXTAR: 
( I2OIOI )]). Because tidal disruption flares last much longer than the 7-ray bursts for which 
Swift was designed, the follow-up (for events involving main sequence stars) can be delayed 
by as much as ~ 1-10 d. 

Applying this analysis to the recently-discovered example of Swift J2058, we find that 
this event is consistent with a tidal disruption of a main sequence star of roughly solar mass 
by a black hole of mass ~ 4 x lO^M©. Its hard spectrum suggests a jet of the sort we have 
described, while its characteristic timescale and X-ray luminosity are likewise consistent 
with this picture. Our model predicts that the jet luminosity should fall when the accretion 
rate begins to fall, but become constant later in the flare when the accretion rate becomes 
sub-Eddington; the observed lightcurve suggests this sort of behavior. Published optical 
photometry is also consistent with our prediction that the optical/UV thermal luminosity 
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should initially be almost independent of time, although this consistency is weakened 
by the large, and plausibly time- dependent, bolometric correction that must be applied. 
Observations using more sensitive detectors might be able to test these predictions more 
definitively. 

Finally, we remark that, although the timescales of Galactic black hole transients are 
such that we can observe a rich phenomenology of spectral states related to the magnitude 
of the accretion rate, the corresponding timescales for supermassive black holes in AGN are 
far too long to permit human observation. Main sequence tidal disruptions, on the other 
hand, present us with a remarkable laboratory for studying how the properties of accretion 
onto supermassive black holes depend on accretion rate. Beginning well above Eddington, 
the accretion rate in these systems declines on timescales from weeks to years, and in a 
fashion we (at least partially) understand. We can hope that further study of these events 
may shed light on the relation between the magnitude of the accretion rate and thermal 
disk emission, jets, and coronal X-rays. 

This work was partially supported by NSF grants AST-0507455 and AST-0908336 
(JHK) and by an ERG advanced research grant and the ISF center for High Energy 
Astrophsics (TP). We thank Michael Stroh for help with the Swift data. 
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